Theory of Rochelle salt: beyond the Mitsui model 



I.V. Stasyuk and O.V. Velychko 

Institute for Condensed Matter Physics 
of National Academy of Sciences of Ukraine, 
1 Svientsitsky Str., 79011 Lviv, Ukraine 

February 2, 2008 



A simple four-sublattice order-disorder model is developed for description of 
phase transitions and dielectric properties of the Rochelle salt crystal. The 
model is developed as a generalization of the semimicroscopic Mitsui model. The 
symmetry properties of lattice and spatial orientations of effective dipoles con- 
nected with the asymmetric structure units in the elementary cell are taken into 
account. The model allows to investigate the temperature and field behaviour 
of transverse (besides longitudinal) components of dielectric susceptibility. The 
infiuence of the transverse electric field E \\ b on the phase transition points and 
spontaneous polarization is studied. 
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1 Dielectric properties of Rochelle salt 

Rochelle salt (RS) is a particular object in the family of ferroelectric crystals with hy- 
drogen bonds. Despite the fact that study of its properties has a long history, the 
structural aspects and mechanisms of phase transitions in this crystal are not conclu- 
sively established. RS becomes ferroelectric (with spontaneous polarization parallel to 
the crystallographic a-axis) in the narrow temperature range between 255 K and 297 K. 
Both nonpolar phases are orthorhombic (P2i2i2i), while the polar phase is monoclinic 
(P2ill). An elementary cell consists of four formula units. 

Numerous data of structural investigations (starting from the early results obtained 
by the X-ray spectroscopy [1] and neutron scattering [2]) do not give a definitive answer to 
the question of microscopic nature of phase transitions in RS. Dielectric relaxation in the 
microwave frequency region and the critical slowing down around the phase transitions 
point to the order- disorder type scenario [3]. Alternatively, the presence of the soft mode, 
which was observed by the far infrared refiectivity and Raman spectroscopy in the lower 
paraelectric phase [4] as well as by microwave dielectric measurements [5] is rather a 
manifestation of the displacive-type transition. 

The soft mode in paraelectric phase is connected with structure changes (such as dis- 
placement of the 0(8) oxygen along the a-axis, rotation of tightly coupled water molecules 
with 0(9) and 0(10) ions) which take place at transition to the ferroelectric phase [6]; it 
is confirmed by the inelastic neutron scattering data [7]. Respective static displacements 
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are the reason of appearance of additional dipole moments of local structure units at 
phase transition to the ferroelectric phase. 

Such displacements can be interpreted also as changes in the population ratio of two 
sites in the disordered paraelectric structure (revealed in the structure investigations 
[8,9]). Large values of the anisotropic temperature factors can be also connected with 
local disorder [10]. The existence of the double atomic positions was taken into account 
in the so-called split-atom model for RS [11]. 

The order- disorder picture of phase transitions in RS forms the basis of the semimicro- 
scopic Mitsui model [12], where asymmetry of occupancy of double local atomic positions 
as well as compensation of electric dipole moments induced in paraphase were taken into 
account. In spite of simplicity (consideration was restricted to two sublattices and in- 
duced local dipole moments were described by means of pseudospins, = ±1/2), the 
model explains quite successfully appearance of two Curie points and effect of deutera- 
tion [12,13]. 

In [14] the model was extended due to inclusion of piezoelectric coupling to the ex- 
ternal field. One should also mention the phenomenological Landau theory [15], adapted 
for systems with a double critical point, which is applicable to the RS crystal in a broad 
range of pressure, substitution concentration of ammonium and temperature. 

The Mitsui model simplifies real structure of the crystal a priori choosing the fer- 
roelectric axis among three two-fold axes thus making an approach essentially "one- 
dimensional" . It is obviously insufficient for more complete description of dielectric prop- 
erties of the RS crystal. We can carry out a generalization, making the model "three- 
dimensional" and taking into account the presence of four (rather than two) translation- 
ally nonequivalent groups of atoms in the unit cell (their positions are mutually connected 
by elements of the point group of the crystal in paraphase [1,2]). Such structure units 
are noncentrosymmetric. Effective dipole moments jli {i = 1, ... ,4) can be assigned to 
them as a whole; the sum of these moments is equal to zero in paraphase. Changes A/Tj 




Figure 1: Orientantions of dipole moments (ordering structure units), producing resulting 
polarization, in the elementary cell of the RS crystal: a comparison between the classical 
Mitsui picture [13] (left) and the proposed approach (right). In the pseudospin formalism 
(S^) = ^{ria — rih), where na j, is a probability (occupation) of respective orientation. 
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in such dipole moments are responsible for appearance of spontaneous polarization in 
the ferroelectric state. Vectors A/Ij are oriented at the certain angles to crystallographic 
axes and possess both longitudinal and transverse components with respect to the a-axis 
(Fig.d). 

Let us use the order-disorder picture for description of such changes. Taking into 
account double equilibrium positions of atoms we come to the effective four-sublattice 
pseudospin model. The model allows to calculate dielectric characteristics in any direction 
and to consider also the effects caused by the influence of the transverse electric field 
(applied perpendicularly to the ferroelectric a-axis). 

In the next section we propose a Hamiltonian obeying symmetry properties of the 
crystal and derive expressions for main thermodynamic characteristics of RS in the mean 
field approximation. The obtained results of consideration of the transverse field effect 
on the polarization and susceptibility are presented in the third section. 
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Figure 2: Orientations of vectors A/ij^ (/c = 1, . . . , 4) in a unit cell of the RS crystal and 
some possible dipole orderings (projection on the plane XY\. (a) high symmetry phase 
(paraphase) - absolute values of pseudospins are equal in all sublattices; (b) ferroelectric 
phase with Ps II -'^ ^ pseudospin values in sublattices 1 and 2 are larger; (c) effect of the 
transverse field - pseudospin values in sublattices 1 and 3 are larger; (d) Ps || X in the 
transverse field - all values are different, sublattices 1 and 2 still prevail. 



3 



2 Four- sublatt ice model: Hamiltonian and thermo- 
dynamics 

According to the above given arguments we take the four-sublattice model as a base 
for simphfied description of phase transitions and dielectric properties of the RS crystal. 
Pseudospin variables S^^, . . . , S^^ describe the before-mentioned changes due to reordering 
of dipole moments of structure units: A/ij^ = d^S^j^. Mean values {S^) = ^{ria — n^) are 
related to differences in position occupancies in the two-minima representation of vectors 
^flik (Fig.El). 

We write down the Hamiltonian of the model in the pseudospin representation: 

ij^j k i,j kj^l 

~ ^ + — S-^ — Si^) — dxEx ^ S-^ 

i ik 

— dyEy — — S-^ + S-^) — d^Ez ''^{Sli — + S^.^ — S-^), (1) 

i i 

where Jkk{hj) and Kki{i,j) describe the inter- and intrasublattice interactions respec- 
tively. The internal field A represents asymmetry in occupancies of the double positions. 
The last three terms in Hamiltonian (P) describe an interaction with components Ea 
{a = X, y, z) of the external electric field. For the sake of simplicity we do not include in 

a term describing tunnelling-like hoppings between equilibrium positions. 

Formula can be considered as a generalization of the Mitsui model Hamiltonian 
[12]: the first four terms are similar to their analogs in that model. Besides the parameter 
r^i, describing the ferroelectric ordering along the a-axis, and the parameter ^, responsible 
for the out of phase ordering of the separated structure elements, there are new parameters 
ri2 and rj^ related to dipole ordering along the h- and c-axes, respectively: 

r^i = \{{Sl) + {8^2) + (SI) + (S^)), e = UiS!) + (^2^) - (SI) - (^D), 

V2 = im) - (SI) - (SI) + (SD), V3 = im) - (SI) + (SI) - (SD). (2) 

Considering Hamiltonian in the mean field approximation we obtain the following 
equations for average values of pseudospins 

(5|) = itanh(i/5i7fe), fc = l,...,4. (3) 

Self-consisted internal fields Hk are given by the expressions 

Hi^2 = {h^ + Ivi) + {h - ^ai^) ± {hy - ^a2r]2) ± {hz + lasr]^), 

H3,4 = iK + ^'ni)-ih-^ai^)Tihy-la2r]2)±ihz + ^a3r]3). (4) 

Here dimensionless quantities h = A/S, ha = daEa/S, 6 = ksT/S, [3 = 1/6, 

ai = [(7^13 + Ku) - ( J + Ku)]/S, a2 = [{Kn - Ku) -{J- Ku)]/S, 

a, = [{Krs - Ku) + (J - Ku)]/S, S = (K,, + Ku) + (J + Ku) (5) 

and the symmetry properties of interaction constants are used. The order parameters 
Vij V2, Vs a^^d the parameter ^ are determined from the set of equations (j2jBl)- Thermo- 
dynamically stable solutions are those with the minimum values of the free energy. In 
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absence of the external field the solution rji ^ 0, ^ 0, r]2 = Vs = ^ corresponds to the 
ferroelectric phase in RS. In this case (S^) = (5*2 ), (5*1) = (S^) and the four-sublattice 
model can be reduced to the Mitsui model. After replacements |?7i i]', |^ —>■ ^' and 
F/ {2N) — > F' /N' in equations (j2}El), one can obtain exactly the same formulae as in that 
case. Nonzero values 172 7^ or 773 7^ are induced in paraphase by the corresponding 
components of the external field. In the ferroelectric phase the parameters 772 and rj^ are 
mutually connected. If one applies electric field along the 6-axis {hy 7^ 0, 772 7^ 0) to 
the RS crystal in the ferroelectric state {rji ^ 0), the third parameter 773 automatically 
becomes nonzero. 



3 Transverse field infiuence on polarization and 
susceptibility 

The main advantage of our model is a possibility to describe the dielectric properties and 
polarization of the RS crystal both along and perpendicularly to the ferroelectric axis. 
There is also an opportunity to consider the effects induced by the external transverse 
field. Fig. |21 illustrates the possible dipole orderings in some important cases when field 
and polarization are parallel to the plane XY {ah). 

Components of tensor of the dielectric susceptibility Xxx = {'2dx/eQVc){drji/dEx) and 
Xyy = {'^dy I eQVc){dri2 / dEy) (where Vc is the unit cell volume) are determined from the set 
of equations (I2H1I) by means of implicit differentiation. Defining Xai3 = i'^dadf^/ SeoVc)Xai3, 
we obtain e.g. for paraphase in presence of the field Ey-. 

_ -Ri(86 — + ^a^Rl 

~ (8e-ii?i)(8e-ia3i?i)-ia3i?r 
_ i?i(86 + ifli-Ri) — iflii?! 

^™ ~ (86 + iai/?i)(8e + ^a2Ri) - \a1a2Rl' ^ ' 

Here Ri = 4(1 — rjl — C,"^) and R2 = —8772^- In the case of the ferroelectric phase the 
parameters r/i and 773 appear besides 772 and ^ in the expressions for Xxx and Xyy Their 
temperature and field dependences are determined from equations (j2Hll)- 

As it follows from the analysis of behaviour of the free energy, the second order type 
of the phase transitions remains unchanged at Ey 7^ 0. In such a case the transition 
temperatures can be determined from the equation 

86 - ^Ri^ (^86 - ia3i?i^ - ^a^Rj = 0, (8) 

which must be solved together with equations following from Q when 771 — 0, 7/3 ~ 
hy7]i 0. 

At small values of the transverse field Ri = 4(1 — ,^q) + o[i?^], R2 ~ Ey, where ^0 
is a solution of the equation ^0 = tanh (/i — |ai^o)]- The critical temperatures shift 
under field proportionally to Ey. Signs and absolute values of the shifts ATc depend on 
magnitudes and signs of the interaction parameters 02 and 03 as well as on the relation 
between them. From pure geometric arguments the parameters Ku and J include in- 
teractions between nearest and next nearest neighbours, respectively. So one can expect 
that K12 > J] it results in the inequality 02 > 03. When we use values Oi = 0.284 and 
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Figure 3: Temperature dependences of average pseudospin values (a), the order parameter 
r]x (proportional to the spontaneous polarization) and longitudinal (c) and transverse (d) 
components of the dielectric susceptibility for different transverse fields at the following 
values of parameters: a\ = 0.284, 02 = 0.1, = —0.25, h = 0.32. 

h = 0.32 chosen to obtain the best fit for critical temperatures of the RS crystal at zero 
field (in this case 5* = 2280 K), the numerical analysis shows that at a-^ < —0.25 the 
ferroelectric region narrows under the field Ey. The direct experimental verification is 
absent, but as some evidence of such a possibility we can consider the results obtained 
in [16,17] by investigations of relaxation phenomena in RS under the external transverse 
field. 

For illustration we give below numerical results for the parameter values a2 = 0.1 and 
03 = —0.25. Temperature dependences of average pseudospin values in the cases depicted 
in Fig. |2l^c,d) are presented in Fig. Efa). Pairs {Sf), (5*1) and (5'|), (S^) demonstrate 
a typical "Mitsui-like" behaviour but the transverse field splits their values creating 
difference between pairs of sublattices even in paraphase. Temperature dependences of 
the order parameter rji, describing spontaneous polarization in the RS crystal, are shown 
in Fig. Efb) for different values of the transverse field. One can see that such a field 
not only narrows the temperature range of the ferroelectric phase but can also suppress 
spontaneous polarization. 

Temperature behaviour of components of the dielectric susceptibility is shown in 
Fig. intc,d). The inverse susceptibility Xxx go^s to zero in the phase transition points 
both at Ey = and Ey ^ (Fig. El^c)). This fact confirms that the phase transitions are 
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of the second order. The transverse component Xyy has jumps in the transition points 
at Ey 7^ 0. Their values are proportional to the second power of the field magnitude 
(Fig. ISfd)). These jumps closely resemble behaviour of the transverse susceptibility of 
the glycinium phosphite crystals in the transverse field [18]. 

Let us make some numerical estimates taking into account the obtained results and 
using the experimental data for Ea and Sb components and for Pg at Ey = 0. The dipole 
moment component can be determined using the maximal value of in the ferro- 
electric phase (-Psimax = 0.25 x 10~^ C/m^ [19])- From the relation Pg = {'^dx/vc)r]i, at 
'7i|max = 0.14 and Vc = 1.04 x 10~^^ cm~^, we obtain dx = 9.26 x 10"^° Cm. Respectively, 
for susceptibility along the a-axis we have Xxx = O-Q^Xxx-, and at Xxx Imax = 3.7 X 10-^ 
it results in Xiximin — 160 (such a value lies inside the experimentally observed range of 
Xxx for ferroelectric phase, see review in [19,20]). 

Estimate for dy component can be obtained using the relation Syy = l+{2dy/ SeoVc)Xyy 
In the ferroelectric phase region Syy = e^, ~ 10 (see [20]; the old experimental data show 
the smooth temperature dependence of Eb in this region). At Xyy = 0.552 we have 
dy = 17.3 X lO"^'' Cm. So the F-component of a dipole moment, connected with reorder- 
ing, is nearly twice as large as the one along the ferroelectric X-axis. 

Let us notice that in this case the field Ey = 18 MV/m corresponds to the value 
hy = 0.01; the shift of T^i is ATd ~ 0.06 K at that field. It means that at the 
fields Ey ^ 1 MV/m the effect will be practically undetectable. The relative change 
of the susceptibility Xyy with temperature in the ferroelectric phase region is also small 
(~ 0.5%). However, the results of numerical estimates can change at another choice 
of the parameter values 02 and 03, so the field effect can be much stronger. As a cer- 
tain argument which points to such a possibility, we can consider the fact, that in the 
GPI crystal (where the effect is caused, as in RS, by the zig-zag geometry of the local 
dipole moment arrangement) the change of Tc is ATc ~ 0.05 K at the transverse field 

~ 1 MV/m. It is obvious that one can check the possibility of such a noticeable effect 
in RS only by direct experimental investigations. 

4 Conclusions 

For description of phase transitions and dielectric properties of the Rochelle salt crystal 
we propose the four-sublattice pseudospin model developed as a generalization of the 
well-known Mitsui model. The introduced model takes into account spatial orientations 
of effective dipoles, which are related to the atomic groups in the unit cell and are 
responsible for the spontaneous polarization. 

The model allows to investigate the temperature behaviour of both longitudinal and 
transverse components of dielectric susceptibility as well as to consider the influence of 
the transverse electric field -E || At the certain relations between the model parameter 
values the increase of this field can lead to the approaching of the lower and higher Curie 
points one to another. Our theory predicts suppression of the spontaneous polarization 
Ps under the field Ey. The effect is similar to the one observed in [16,17]. The changes 
in the transverse susceptibility Xyy the transition points are similar to the phenom- 
ena detected in the GPI crystal [18]. As it follows from our consideration, there should 
exist jumps in Xyy which increase proportionally to Ey. However, performed here nu- 
merical estimates indicate that these effects could have small magnitudes. Only future 
experimental studies can give a final answer. 
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The ideas forming the basis of the four-sublattice model can be used for description of 
the mixed system RSi_a;-ARSj;, where the significant changes in temperature behaviour 
of the longitudinal and transverse dielectric susceptibilities are observed at increase of the 
concentration x of ammonia groups and in the high concentration region (0.89 < x <1) 
a polar phase appears with polarization along the 6-axis [21,22]. 
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